Dynamics of the superconducting condensate in the presence of a magnetic field. 
Channelling of vortices in superconducting strips at high currents 
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On the basis of the time-dependent Ginzburg-Landau equation we studied the dynamics of the 
superconducting condensate in a wide two-dimensional sample in the presence of a perpendicular 
magnetic field and applied current. We could identify two critical currents: the current at which the 
pure superconducting state becomes unstable (J c 2^) and the current at which the system transits 
from the resistive state to the superconducting state (J c i < J C 2)- The current J C 2 decreases mono- 
tonically with external magnetic field, while J c i exhibits a maximum at H* . For sufficient large 
magnetic fields the hysteresis disappears and J c i = Jci = Jc- In this high magnetic field region 
and for currents close to J c the voltage appears as a result of the motion of separate vortices. With 
increasing current the moving vortices form 'channels' with suppressed order parameter along which 
the vortices can move very fast. This leads to a sharp increase of the voltage. These 'channels' 
resemble in some respect the phase slip lines which occur at zero magnetic field. 
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It is well-known that the resistive state in supercon- 
ducting wires or stripes with diameter/ width less or com- 
parable to the coherence length £ is realized through the 
appearance of phase slip centers2^i4. The phenomenolog- 
ical theory of phase slip centers (PSC) was first proposed 
in Ref A. According to this theory a PSC is a region with 
size of order £ where the order parameter is strongly sup- 
pressed. The normal current density produced by the os- 
cillation of the order parameter in the phase slip center 
decays on a larger distance scale Aq ^> 

Current- voltage characteristics of such a system is usu- 
ally irreversible^. It is possible to distinguish two crit- 
ical current densities: j c i - current density at which the 
superconducting state becomes unstable and j c \ < jc2 - 
current density below which the phase slip solution does 
not exist in the system^. In recent work£ it was claimed 
that the current j c ± is defined by the competition between 
two relaxation times: the relaxation time of the absolute 
value of the order parameter tui and the relaxation time 
of the phase of the order parameter r^. The phase slip 
solution does exist when roughly t\^\ > t$. 

Because an applied magnetic field suppresses the or- 
der parameter and leads to the appearance of screening 
currents in the sample it is naturally to expect that it 
affects the phase slip process in the superconductor. For 
example in Refs. 6 ' 7 it was shown that a parallel (to the 
direction of the injected current) magnetic field modifies 
the critical currents j c ±, j C 2» and the stair structure of 
the current-voltage characteristic 6 . 

If we apply a perpendicular magnetic field the situa- 
tion becomes more complicated. In this case the screen- 
ing currents induced by the magnetic field decreases the 
current density j on one side of the superconductor and 
increases j on the other side of the sample (see Fig. 
1) and the current density in some part of the stripe 
becomes smaller than j c i (if total current is equal to 
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FIG. 1: Schematic current density distribution in a supercon- 
ducting stripe in a perpendicular magnetic field at the critical 
currents I c i(H) (solid lines) and I C 2(H) (dotted lines). In the 
insert a schematic view of the considered set up is shown. 



J = jci IF). In accordance with Refi, in that part of 
the sample the phase slip process cannot be realized. We 
should increase the applied current in order to satisfy 
the condition j > j c \ in any point of the line connecting 
the two opposite sides of the strip. Further increasing 
H the slope of j(x) increases (at current J = J c i(H)) 
and at the moment when the current density on the left 
side of the sample reaches j C 2 (see Fig. 1) both critical 
currents becomes equal to each other J c2 = J c \ at field 
H = H*~(j c2 -j cl )/W. 

To check the above predictions we studied the 
current-voltage characteristics of quasi-two-dimensional 
superconductors using the generalized time-dependent 
Ginzburg-Landau (TDGL) equation. The latter was first 
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FIG. 2: Current-voltage characteristics of a superconducting 
strip in a perpendicular magnetic field for sweep up and sweep 
down regimes. 
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FIG. 3: Dependence of the critical currents j c i and j C 2 on the 
applied perpendicular magnetic field. 



written down in the work of Ref, 8 

= (V-iA)V + (l-|^| 2 )V'- (1) 

and should be supplemented with the equation for the 
electrostatic potential 

Atp = div(Irru>*(V-iA)V>)), (2) 

which is nothing else than the condition for the conser- 
vation of the total current in the wire, i.e. divj = 0. 
In Eqs. (1,2) all the physical quantities (order parameter 
ijj = \ip\e 1 ^, vector potential A and electrostatic potential 
tp) are measured in dimensionless units: the vector po- 
tential A and momentum of superconducting condensate 
p = V0 — A is scaled by $o/(27r£) (where $o is the quan- 
tum of magnetic flux), the order parameter is in units 
of Ao and the coordinates are in units of the coherence 
length £(T). In these units the magnetic field is scaled 
by H C 2 and the current density by jo — c^q/8tt 2 A 2 ^. 
Time is in units of the Ginzburg-Landau relaxation time 
tql = 47rcr Il A 2 /c 2 = 2Th/nAl, and the electrostatic po- 
tential {ip) is in units of ipo = c$o/8ir 2 t;\(7 n = h/2eTGL 
(a n is the normal-state conductivity). The parame- 
ter u is about 5.79 according to Ref£. We also put 
A = (0, Hx,0) in Eq. (1,2) because we considered the 
case when the effect of the current-induced magnetic field 
is negligible. 

In our calculations we mainly used vacuum- 
superconductor boundary conditions (V — iA)iJj\ n = 
and V<p|„ = except for the regions where current was 
injected (see insert in Fig. 1). In those points we used 
the normal metal-superconductor boundary conditions 
tp = and Vtp\ n = —j. 

We found two jumps in the current-voltage character- 
istics (see Fig. 2). The first one is connected with the 
transition of the sample from the supercondicting to the 



resistive state (but with \ip\ ^ 0). The second transition 
is the one from the resistive to the normal state. The 
hysteresis connected to the second transition survives till 
very high magnetic field. The latter is connected with the 
stability of the normal state carrying current at tempera- 
tures lower than the critical temperature (see for example 
Ref£). For such a system it was shown that the normal 
state may exist (at T < T c ) till very low current densi- 
ties if there is no finite superconducting nucleus. In our 
case we inject current through part of the cross-section 
of the sample (see inset of Fig. 1) and in this way we 
provide a more optimal condition for the nucleation of 
superconductivity in the corners of the superconducting 
sample where the current density is minimal even in nor- 
mal state. We want to stress that the second hysteresis 
crucially depends on the geometrical parameters of the 
sample and it is not the subject of the present paper to 
study that effect. We only study here the first hysteresis 
which is almost sample-independent. 

In Fig. 3 we present the dependencies of J c i(H) and 
Jc2{H) for different widths of the sample. We found that 
in accordance with our predictions the current J c \ in- 
creases at low magnetic fields and at some H = H* it 
becomes equal to J c2 . Beyond H* that type of hysteresis 
no longer exists in the system and the critical current 
J c = J c i = J C 2 decreases with increasing magnetic field. 
Only for relatively narrow samples with W ~ £ there is 
a second peak in the J C {H) dependence. The reason for 
that peak is probably connected with strong nonlinear 
effects. Indeed, the current density is related to the mo- 
mentum of the superconducting condensate (p = V0— A) 
as j = p(l — p 2 ) in accordance to the Ginzburg-Landau 
relation. The second peak occurs approximately at a 
field H s when the first vortices penetrate the sample in 
absence of injected current. At H = H s the value of p 
on the edges is close to unity for samples with W ~ £. 
It means that the term (1 — p 2 ) may be very important 
in that range of fields. In wider stripes the above effect 
is negligible and there is no second peak in the J C (H) 
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curve. 

For magnetic fields larger than H s , vortices will pen- 
etrate the superconductor and occupy the central part 
of the sampl oV . If one injects current into the sam- 
ple the vortex-filled region starts to move to one side of 
the sample and when the vortex dome touches the sam- 
ple boundary the vortices will leave the sampl o?i 10 . This 
value of the current is the critical one. It is essential 
that the distribution of the current density in the vortex 
dome is practically constant in the resistive state (in nar- 
row samples with W < A = \ 2 /d) even in the absence 
of bulk pinningii. At fields H ;§> H s the vortex dome 
occupies almost the whole sample and it means that at 
J ~ Jd(0) = jdW the conditions for the occurrence of 
phase slip lines (or vortex channels) will be fulfilled and 
it will lead to the appearance of a stair structure in the 
I-V characteristics and to a sharp increase of the voltage 
at J > J cl (0)A£. It also follows that at relatively large 
H this current should be field- independent. The lower 
boundary of the field-independent region is defined by 
the condition that the current density distribution should 
be almost uniform (Hi ower ~ H s ). The upper boundary 
H U pp e r depends on how large is magnetic field and how 
strong the order parameter is suppressed in vortex filled 
region because according to Refsi&£ both these factors 
affect the value of j c \ (it leads to a decrease of it). 

We found that the voltage sharply increases at J > J c 
when vortex slip lines or vortex channels appear in the 
sample (compare Fig. 2 and Fig. 4 at H=0.3 and H=0.6). 
In contradiction with the above prediction the current 
at which the vortex channelling occurs decreases with 
increasing H (see Fig. 2). The most probable reason is 
that the current density strongly varies over the width 
of the sample for H ;§> H s in case W = 9£ because even 
close to H C 2 there will be only three rows of vortices (see 
Fig. 4). 

Our results coincide qualitatively with recent experi- 
mental results on the I-V characteristics of carbon nan- 
otube (compare Fig. 13 in Ref. 13 with Fig. 3 in our pa- 
per) measured at different magnetic fields. In that work 
a non-monotonous dependence of J c \ on H was observed 
at low magnetic fields. In Refii^ a stair structure in I- 
V characteristics of Mo/Si multi-layer stripes was found 
which appears practically at the same value of the in- 
jected current, as predicted in this paper, in wide field 
region. 

In conclusion, we found that an external magnetic field 
strongly affects the resistive state of mesoscopic wires and 



stripes. It leads to a shrinking of the hysteresis in the 
current-voltage characteristic at relatively high values of 
H. If the magnetic field is perpendicular to the current di- 
rection the dependence J c \ (H) is a non-monotonous func- 
tion of H at low magnetic field. At high magnetic fields 
J c i and J C 2 monotonically decreases with increasing H 
for both perpendicular and parallel^ orientation of the 
applied magnetic field. 
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FIG. 4: Contour plot of the order parameter in the resistive 
state of a superconducting strip at different values of the mag- 
netic field and applied current. The width of the sample is 
9£, the length is 40 £, and we used the parameter 7 = 10, 

(i> = J/UoW). 



* Electronic address: denis@uia.ua.ac.be 
t Electronic address: peeters@uia.ua.ac.be 

1 Everywhere in the paper we define the current by J = 
I/d = f Q W j(x)dx. 

2 B.I. Ivlev and N.B. Kopnin, Adv. Phys. 33 (1984) 80. 

3 M. Thinkham, Introduction to superconductivity, 
(McGraw-Hill, NY, 1996), Chpt. 11. 



4 R. Tidecks, Current-induced nonequilibrium phenomena in 
quasi-one-dimensional superconductors, (Springer- Verlag, 
Berlin, 1990). 

5 W. J. Skocpol, M. R. Beasley, and M. Tinkham, J. Low 
Temp. Phys. 16 (1974) 145. 

6 A.M. Kadin, W.J. Skocpol, and M. Tinkham, J. Low 
Temp. Phys. 33 (1978) 33. 



4 



7 S. Michotte, S. Matefi- Tempfli, L. Pi raux, D. Y. Vodola- 
zov, and F. M. Peeters, cond-mat/0309699 

8 L. Kramer and R.J. Watts-Tobin, Phys. Rev. Lett. 40 
(1978) 1041. 

9 M.Yu. Kypriyanov and K.K. Likharev, Sov. Phys. Solid 
State 16 (1975) 1835. 

M. Benkraouda and J.R. Clem, Phys. Rev. B 58 (1998) 
15103. 

1 G. M. Maksimova, N. V. Zhelezina and I. L. Maksimov, 
Europhys. Lett. 53 (2001) 639. 

2 The coexistance of phase slip process and flux flow was 
discussed already in the paper: S.V. Lempitsky, Zh. Eksp. 
Teor. Fiz. 90 (1986) 793. The author predicted (using 
the modified London approach) that the resistive state 



in wide (W S> £) stripes may be realized at nonzero 
magnetic field for low currents through flux flow and for 
high currents through the appearence of phase slip lines. 
The nonmonotonous dependence of J c i(H) and the field- 
independent current at which the phase slip lines (vortex 
channels) appear at high H were not discussed. 

A. Kasumov, M. Kocial, M.Ferrier, R. Deblock, S. Gueron, 

B. Reulet, I. Khodos, O.Stephan, and H. Bouchiat, 
preprint. 

A.Y. Aladyshkin, A.Y. Klimov, M.A. Kucherenko, V.V. 
Kurin, V.L. Mironov, Y.N. Nozdrin, and E.E. Pestov, un- 
published. 



